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Velocity Distributions of Tracer Particles in Thermal Counterflow in Superfluid 4 He 
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Quantum turbulence accompanying thermal counterflow in superfluid 4 He was recently visualized 
by the Maryland group, using micron-sized tracer particles of solid hydrogen (J. Phys. Soc. Jpn. 
77, 111007 (2008)) . In order to understand the observations we formulate the coupled dynamics of 
fine particles and quantized vortices, in the presence of a relative motion of the normal and super- 
fluid components. Numerical simulations based on this formulation are shown to agree reasonably 
well with experimental observations of the velocity distributions of the tracer particles in thermal 
counterflow. 
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PACS numbers: 67.25.dk, 67.30.eh 



I. INTRODUCTION 



Quantum turbulence is a type of turbulence that is 
strongly affected by quantum effects. It is observed in 
superfluids, in which rotational flow of the superfluid 
component can take place only through the presence of 
quantized vortex filaments. Turbulence in the superfluid 
component then takes the form of a tangle of these quan- 
tized vortices^. Quantum turbulence was first observed 
and studied 3 - in thermal counterflow in superfluid 4 He. 
Such counterflow is associated with two-fluid behaviour, 
in which the superfluid component, density p s , coexists 
with a normal component, density p n , which carries all 
the entropy in the system 4 -. The two components can 
have separate velocity fields, v s and v n . Heat is carried 
in such a system by counterflow of the two fluids, with no 
net mass flow, the superfluid component moving towards 
the source of heat and the normal component away from 
it. Above a certain critical heat current this counterflow 
proves to be unstable, the superfluid component becom- 
ing turbulent. The normal fluid may also become turbu- 
lent, but probably only at heat currents larger than those 
with which we are concerned in this paper. This type of 
counterflow turbulence has no classical analogue. The 
vortex lines present in the turbulent superfluid compo- 
nent interact with the normal fluid, giving rise the force 
of "mutual friction" between the two fluids, and study of 
the behaviour of this mutual friction provided the first 
experimental evidence for the existence of counterflow 
turbulence^. Since then, many experimental, theoreti- 
cal and numerical studies of counterflow turbulence have 
been reported. Schwarz carried out numerical simula- 
tions with the vortex filament model to confirm the ex- 
istence of a self-sustained statistically steady turbulent 
state in the superfluid component^, and a more satisfac- 
tory simulation based on a full Biot-Savart simulation 



has recently been published by Adachi et al~. However, 
until recently there has been no direct visualization of 
the flow, such as might convince a skeptic that thermal 
counterflow is a reality and might provide detailed in- 
formation about the states of turbulent flow. We shall 
first describe briefly the two visualization techniques that 
have been developed recently, and then focus on the in- 
terpretation of one particular set of experimental results. 

The first technique relies on the production within the 
helium of micron-sized tracer particles formed from solid 
hydrogen, and motion of these particles has been fol- 
lowed by both Particle Image Velocimetryl and Particle 
Tracking Velocimetry^. The second technique^ uses 
metastable triplet state He2 molecules as tracers. These 
molecules are not trapped by vortex lines above IK, and 
they respond only to motion of the normal fluid; they 
have been used recently to show that at a sufficiently 
large heat flux the normal fluid in thermal counterflow 
does become turbulent. In this paper we shall be con- 
cerned with the interpretation of experiments on ther- 
mal counterflow that use as tracers particles of hydrogen; 
these tracers can be trapped (intermittently) by vortex 
lines above IK, so that interpretation is less straightfor- 
ward than is the case for the molecules. However, the 
experiments by Paoletti et al£ with which we shall be 
concerned relate to heat fluxes that are rather small, so 
that we are at least spared the need to consider the effect 
of a turbulent normal fluid. 

In the experiments by Paoletti et al£ the tracer par- 
ticles in a thermal counterflow were observed to divide 
into two groups: those that move freely in the direction 
of the normal fluid, and those that are trapped into vor- 
tex lines and move in the opposite direction. The velocity 
distribution of the trapped particles is broader, reflecting 
an irregular motion of the vortices. This difference was 
visible in the observed trajectories: smooth for the first 
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group but more irregular for the second. The primary 
aim of this paper is to account for Paoletti's observations 
by modelling the coupled dynamics of particles and vor- 
tices and using the model in numerical simulations. The 
motion of a free particle that is not trapped by vortices 
is determined by a Stokes drag from the viscous normal 
fluid and by inertial effects which arise from both super- 
fluid and normal fluid. The motion of a particle trapped 
by a vortex is determined by the Stokes drag and iner- 
tial effects, together with the effect of the tension in the 
vortex to which the particle is attached. The motion of a 
vortex line depends on its shape and on the force of mu- 
tual friction with the normal fluid, but it can be modified 
by the presence of the particle. The problem of particle- 
vortex interaction can be tackled at different levels of 
approximation. At the most microscopic leveliir— , the 
close interaction of particle and vortex and the trapping 
of the particle onto the vortex require the calculation 
of a complex system of image vortices (mathematically, 
trapping is the reconnection of the vortex with its own 
image inside the particle). Furthermore, the particles of 
hydrogen are of irregular and unknown shape, so that 
this kind of calculation is hardly relevant in practice. As 
described in Section II, we have therefore used a simple 
model which neglects the details of the particle-vortex 
interaction that occur on length scales of the order of 
the particle size or less, and which is sufficiently simple 
to allow the study at a useful and instructive level of the 
coupled dynamics of many particles and many vortices si- 
multaneously. In Section III we present an even simpler 
model, applicable only at the smallest velocities, when 
the configuration of vortex line is modified to a negligi- 
ble extent by the presence of the particles. We summarise 
our results in Section IV. 



II. COMPUTATIONS 

We use the vortex filament model£&i£ and represent 
the vortex as strings of discrete points s (£,£), where t 
is time and £ is arc length. The superfluid velocity at 
a point r due to a filament is given by the Biot-Savart 
expression: 

K f (si - r) x dsi 
v s (r) = — ■ -3 — , 1 

where k is the quantum of circulation, Si is a point on 
the filament, and integration is taken along the filament. 
Calculating the velocity v s at a point r = s on the fila- 
ment causes the integral to diverge as S\ — > s. To avoid 
this divergence, we separate the velocity s of the filament 
at the point s into two components^: 

k , /2(M ) 1/2 \ k ['(s-r)xds 

s = — s' x s" In v T,„ ; H / — . 

4tt ^ e V2£ J AttJc \s - r\ 3 

(2) 

The first term is the localized induction field arising from 
a curved line element acting on itself, and Z+ and I- are 



the lengths of the two adjacent line elements after dis- 
cretization of A£, separated by the point s. The prime 
denotes differentiation with respect to the arc length £. 
s', s" , and s' x s" are the tangent vector, the normal 
vector, and binormal vector, respectively. The parame- 
ter £0 is the cutoff corresponding to the core radius. The 
second term represents the non-local field obtained by in- 
tegrating Eq. (|T|) along the rest of the filament, excluding 
the neighborhood of s. We take account of the contribu- 
tions from the boundary-induced velocity field v St ^ and 
applied field v s ^ a in Eq. @. 
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FIG. 1: (Color online) (a) Illustrating the way in which a 
vortex, moving relative to the normal fluid, is distorted by 
a trapped particle; a vortex ring, instantaneously in the x- 
y plane, propagates along the z-axis, with and without a 
trapped particle, the normal fluid being at rest ; the ring 
gradually collapses, under the combined influence of mutual 
friction and the viscous drag on the particle. (The (red) cusp 
to which the particle is attached points in the direction of 
—z.) (b) and (c) Schematic of the tension force acting on an 
element of vortex line, with and without an attached particle. 
We emphasize that (c) is schematic only; the real vortex con- 
figuration in (c) requires a representation in three dimensions, 
as in (a). 

We now explain how we handle to motion of a vor- 
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tex line to which a particle is attached at a particular 
point. The way in which a vortex line is distorted by a 
trapped particle is illustrated in Fig. 1(a). We shall as- 
sume that the vortex line is smooth on a scale of order the 
separation between the discrete points used in the com- 
putation, except in the immediate neighbourhood of the 
particle, and that the particle is situated at one of these 
points s(£) (Fig. 1(c)). The forces acting on the particle 
are: a net force due to the tension in the sections of vor- 
tex line that adjoin the particle, which is related to the 
angle between these two sections of line; the drag force 
on the particle due to the viscosity of the normal fluid 
flowing relative to the particle; and inertial forces aris- 
ing from both the superfluid component and the normal 
fluid; and a Magnus force. We note that the angle be- 
tween the sections of vortex line adjoining the particle is 
not equal to the angle between the lines joining the parti- 
cle point to the adjacent discrete vortex points, s(£ + A£) 
and s(£ — A£). Furthermore, we do not know the precise 
shape of the vortex in the immediate vicinity of the parti- 
cle. This potential difficulty can be handled in a number 
of different ways. Here we shall handle it by consider- 
ing the forces acting on the whole section of vortex line 
extending between the two mid points s(£ + A£/2) and 
s(£ — A£/2), including the particle, regarded as a single 
entity. 

To make our procedure more clear we shall first con- 
sider the forces acting on this section of vortex line in 
the absence of the particle; in this case there can be no 
kink in the vortex line at s(£) (Fig. 1(b)). First there is a 
force, Ft due to the tension in the line, given by 



F T = T [s'(£ + A£/2) -a'(£- A£/2)], 
where To is the line tension, given by 

'2(/ + /_) 1/2 \ 



(3) 



(4) 



not yield any contribution from the non-local field to the 
component of s in a direction along the vortex line; this 
component does not have any real physical effect on the 
vortex dynamics, although it can lead to a motion of the 
discrete vortex points. 

Our analysis so far has assumed that the temperature 
is vanishingly small. At a finite temperature a force of 
mutual friction, given by 



MF 



(-70s' x [*' x (v n - s)] -j' Q s'x {v n - s)) A£, (7) 



must be added to the balance of forces, where 70 and 7 
are temperature-dependent parameters. Therefore s is 
now given by 



Fm + Fj 



MF 



0. 



(8) 



It is easily shown that s is then given by 



s = s + as' x (v n - so) - as' x \s' x (v n - so)], (9) 

where s is given by Eq.©, and a and a' are the more 
usual mutual-friction parameters. 

Now consider the behaviour of a spherical tracer par- 
ticle, radius a p (< 0.5A£) and volume V p . The density 
of the material of the particle is assumed to be the same 
as that of the helium, so that the particle is neutrally 
buoyant. Consider first the motion of the particle when 
it is not trapped on a vortex. Its equation of motion is 



3 dv p 
2 PV ^ 



fs + fi. 



(10) 



We have allowed for the fact that the effective mass of 
the particle is enhanced in the fluid by an amount equal 
to half the mass of fluid displaced. The force fs is the 
Stokes drag associated with the viscosity of the normal 
fluid, given by 



In a computation we make the approximation that s'(£ + 
A£/2) is equal to the slope of the straight line joining the 
points s(£) and s(£ + A£). Secondly there is a Magnus 
force, given by 



F M = p s ns' x [s(£) - u s ,„i]A£, 



where v s n \ is the second term in Eq.([2]) evaluated at the 
point s(£). Since the mass of the core of a vortex is 
negligibly small, these two forces must balance, provided 
that there is no mutual friction. Therefore the velocity s 
is given by 



M 



0. 



(6) 



If the shape of the vortex is sufficiently smooth, and if 
A£ is sufficiently small, the first term in Eq.© is equal 
to Tos"A£, so that Eq.© yields a velocity s(£) similar to 
that given by Eq. ([2]) , but derived in a different way. The 
only difference is that Eq.©, in contrast to Eq.@, does 



67ra p ^„(t) n - v p ), 



(11) 



where fi n is the viscosity of the normal fluid. We have 
assumed that the particle Reynolds number is small, and 
that the mean free path of the excitations constituting 



(5) the normal fluid is much smaller than a p . The force f\ is 



an inertial force acting on the particle associated with any 
acceleration of the normal and superfluid components, 
and it is given by 



3 Dv s 
fl = 2^ P ~DF 



3 T . Dv n 



(12) 



In our simulations we neglect the inertial force on the 
normal fluid, because v n is assumed to be constant in 
time and spatially uniform. 

Finally we consider the situation when the tracer par- 
ticle is trapped on a vortex at the point s(£) (Fig. 1(c)). 
The element of vortex between the points s(£ + A£/2) 
and s(£ — A£/2) is now subject to the extra forces fs 
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and fi, and the total force on this element must lead to 
the particle acceleration dv p /dt; i.e. 

lpV P ^ = F T + F M + F M F + fs + fi- (13) 

It is important to understand that in calculating the 
forces Ft, Fm and Fmf we must take the element of 
vortex between the points s(£ + A£/2) and s(£ — A£/2) 
to have the shape that it has when the particle is at- 
tached; this shape is not known in detail, and it must 
involve, for example, a discontinuity in the slope s' at 
the position of the particle. To make further progress we 
must therefore make approximations. We shall make two 
assumptions, which cannot be accurately correct; we can 
hope that they do not introduce serious errors, although, 
as we shall see, there is some evidence that they do lead 
to small errors. First, we shall assume that the element 
of vortex moves with an average velocity equal to that of 
the particle, v p . And secondly we shall assume that the 
effective length of the element, relevant to the calculation 
of Fm and Fmf, is not significantly different from that 
shown by the smooth curve in Fig. 1(c). We judge that 
these assumptions are likely to be good if the particle is 
significantly smaller than A£. Then the values of Fm and 
Fmf are given by Eqs.© and J7]), but with s replaced 
by v p , while Ft is still given by Eq.Q. Therefore 

3 m dv P a i \ , 3 t/ DVs 

2 P p ~dt = 6?ra P^ n ( Vn ~ v p) + 2 Ps p_ Di" 
T [ S '(£ + A£/2)- S '(£-A£/2)] + 

P s ks' x (v p - V s , nl )A£ + 

(70s' x [«' x (« p - v n )] + To V x (v p - v n )) AC (14) 

We make the important remark that the acceleration 
given by Eq. (|14p ought not to depend on the magnitude 
of A£. A reduction in A£ reduces the magnitudes of the 
last two terms on the right hand side of this equation 
(the Magnus force and the force of mutual friction), both 
of which are proportional to A£, but this reduction is 
cancelled by an increase in tension term To[s'(£+A£/2) — 
s'^-A^/2)}. 

We note that, correctly, the component of i> Si „i parallel 
to the vortex makes no contribution to the velocity of the 
particle, although, as we have already noted, it does make 
a contribution to the motion of the discrete vortex points 
that are not attached to a particle. 

In summary, then, the dynamics of the untrapped par- 
ticles is described by Eq. fTU)) . and coupled dynamics of 
the vortices and trapped particle is described by Eq. (fT4"| , 
where v Sjn i is given by the second term in Eq.@. 

In practice, the terms involving accelerations in 
Eq. (fl"4f are small and can be neglected, so that the size 
of the particle enters only in the Stokes drag as a p . In 
our simulations we have taken as a p = 1.0 x 1CP 4 cm. 
We have used two values of the numerical space resolu- 
tion A£: 2.0 x 10~ 4 cm and 8.0 x 10~ 4 cm. Our time 
resolution is At — 1.0 x 10~ 6 s. We use periodic bound- 
ary conditions with a box of size 0.1 3 cm 3 . We have 



used three temperatures: 1.9 K, 2.0 K, and 2.1 K. The 
velocity of the normal fluid is assumed constant in time 
and spatially uniform. First we make thermal counter- 
flow quantum turbulence that is statistically steady and 
satisfies the well-known relation L = 7 2 v 1 i S where L is 
the vortex line density and 7 is a temperature-dependent 
parameter—. After that, we place 120 particles trapped 
on the vortices (Fig. [5]). The turbulence continues to 
evolve in a steady state, and we study the distribution of 
particle velocities. A particle velocity is obtained from 
the distance moved by the particle in a fixed time inter- 
val, taken to be 10 ms when A£ = 2.0 x 10~ 4 cm and 100 
ms when A£ = 8.0 x 10~ 4 cm. Some particles become 
detrapped, and these are ignored; this detrapping will be 
considered in a later paper. 

We note that our values of A£ are not large compared 
with a p , contrary to an assumption made in our deriva- 
tion of Eq. (fT4")) . These small values were chosen in order 
to ensure as far as we can that our simulations describe 
correctly the evolution of sharp bends in the vortices. 
Strictly speaking, therefore, our simulations relate to par- 
ticles that are significantly smaller than 10 -4 cm in ra- 
dius, but which experience a Stokes drag equal to that 
experienced by the larger particle. Although smaller par- 
ticles are likely to be detrapped more easily, we can hope 
that the behaviour of the trapped particles is not signif- 
icantly affected. 




FIG. 2: (Color online) A snapshot of the coupled dynamics 
of the tracer particles and the vortices at T=1.9 K. A£ = 
2.0 x 10 -4 cm. All of the 120 particles are trapped by vortices. 
The direction of the normal flow v n = 0.30 cm/s is upward, 
and that of superflow v a given by Eq. (|15[) is downward. The 
vortex line density is L ~ 6000 cm -2 . 



In Fig|3] we show examples of the components of the 
velocity distributions of the tracer particles parallel to 
velocity of the normal fluid. Superfluid velocities are cal- 
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culated from the equation 
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(15) 
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FIG. 4: (Color on line) Illustrating the sliding of a particle 
along a vortex. 



FIG. 3: (Color online) Velocity distributions of the parti- 
cles in the direction of the flow of the the normal fluid. 
A£ = 2.0 x 1CP 4 cm. The velocity v p is normalized by the 
normal fluid velocity v n . The red, green and blue points and 
lines are at T=1.9 K, 2.0 K, 2.1 K, respectively, and the cor- 
responding normal fluid velocities are 0.30, 0.20 and 0.09 cm 
s~ , respectively (corresponding to approximately the same 
vortex line density). The vertical lines refer to the velocity of 
the superfluid given by Eq. (fl5|) . 



We see that the velocity of the particles is significantly 
slower than that of the vortices, which move at a velocity 
that is only slightly smaller than that of the superfluid. 
This feature is probably associated with a sliding of the 
particles along the vortices, as we shall see more clearly 
in the next Section. 



III. BEHAVIOUR OF THE TRAPPED 
PARTICLES IN THE LIMIT OF SMALL 
VELOCITIES 

When velocities are sufficiently small, the Stokes drag 
on a particle is insufficient to cause significant distor- 
tion of the vortices. In that case the velocity with which 
a particle moves can be calculated by assuming that it 
simply slides along the vortex to which it is attached, 
the properties of the vortex array being determined from 
simulations carried out in the absence of particles. 

Consider the behaviour of one particle under these con- 
ditions on a vortex line that is inclined at an angle 6 to 
the direction of flow of the normal fluid (FigBJ . We need 
first to calculate the velocity v p of the particle along the 
vortex, when the vortex is itself moving with the veloc- 
ity «l, given by Eq.Q, evaluated in the absence of any 
trapped particles; «l is not necessarily in the plane de- 
fined by the vectors v p and v n . The Stokes force acting 



on the particle is given by the vector equation 

F = 67r// n a p (v n - v p - v L j_), (16) 

where vlj_ is the component of vl perpendicular to the 
vortex line. The component of this force along the vortex 
line is 

F L = 6-Kfi n a p (v n cos# - Up). (17) 
This component must vanish. Therefore 

Vp — Vn COS0. (18) 

Therefore the component of the particle velocity in the 
direction of v n is 

v p cos 9 + t! L _l_ n = v n cos 2 9 + V L±n , (19) 

where «L± n is the component of vl_l parallel to v n . 

We take an average over all particles, assumed to be 
randomly distributed on the vortex lines. The average 
component of the particle velocity in a direction perpen- 
dicular to v n must vanish, and therefore we can write for 
the average particle velocity, which is in the direction of 

Vn, 



(v p ) = v n (cos 2 9) + (Vh± a ). 



(20) 



The average («L±n) can be evaluated from the results of 
simulations of the counterflow vortex tangle^ 

a_Ln> = -7^— / Vn • [s' X (s' X S )]g^. (21) 



Symmetry requires that the vector ([s' x (s' x s)]) is par- 
allel to v n , so that 

(vLXn) = -^J I 8 ' >< ( S ' >< SM- (22) 
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Evaluation of the average (cos 2 9) requires a knowledge 
of the anisotropy of the vortex tangle. In fact it is related 
the anisotropy parameter, introduced by Schwarz and 
evaluated in a full Biot-Savart simulation by Adachi et 



FIG. 5: (Color on line) The results of our theories and com- 
putations. The solid line shows v p — v s - 



J„ = ^ /(I 



so that 



(23) 



(24) 



Thus we have for the mean particle velocity 



(Vp) =Wn(l-/||) + |(v L ±n>| = f-^(l-/||) + ^^ 

V p n V s 

(25) 

where we have used Eq. (fT5"|) . Eq. ([2"5j) is our hnal result 
for the mean particle velocity in the limit of small ve- 
locities; The integral 7|| has already been calculated in 
the simulations of Adachi et a££, and we have now cal- 
culated the quantity (vLj_ n ) from the results of the same 
simulations. 



IV. SUMMARY OF OUR RESULTS 
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The results of our various simulations and calculations 
are displayed in Fig[5] The points relating to our simula- 
tions (Section II) are the particle velocities at the peaks 
of distributions of the type shown in FigEl and the error 
bars are measures of the standard deviations of these dis- 
tributions. We have chosen to display them in this way 
to facilitate comparison with the experimental results of 
Paoletti et aL— , which are displayed in the same way in 
their Fig. 7(b). We see that within the large scatter, seen 
in both the results of our simulations and the experimen- 
tal results, the different computed velocities and the ex- 
perimentally observed velocities are in rough agreement. 
The velocities are in the direction of, but significantly less 
than, v e , but are probably proportional to v s . We judge 
that this proportionality is a consequence of a simple slid- 
ing of the particles along the vortices, without significant 
distortion of the vortex tangle. If this is indeed the case, 
then we see that our simulations probably do not extend 
to velocities at which there is significant distortion; judg- 
ing from the experimental results, distortion sets in only 
at superfluid velocities above roughly 0.4 cm s _1 . Al- 
though, as we say, our results are in rough agreement 
with experiment, a careful examination of Fig J5] suggests 
that agreement is not as good as we might have hoped. 
Our three different sets of theoretical results (relating to 
two different values of A£, and to the calculation of Sec- 
tion III) yield values of the mean particle velocities that 
seem to differ from one another by small but significant 
amounts (the standard deviations of the means should be 
much less that those shown in FiglSJ). These discrepan- 
cies may well be associated with our failure to take into 
account adequately the precise shape of a vortex in the 
immediate vicinity of a particle. 
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